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ABSTRACT 


A physical inversion scheme has been developed, dealing with cloudy as well as cloud-free 
radiance observed with ultraspectral infrared sounders, to simultaneously retrieve surface, 
atmospheric thermodynamic, and cloud microphysical parameters. A fast radiative transfer 
model, which applies to the clouded atmosphere, is used for atmospheric profile and cloud 
parameter retrieval. A one-dimensional (1-d) variational multi-variable inversion solution is 
used to improve an iterative background state defined by an eigenvector-regression-retrieval. 
The solution is iterated in order to account for non-linearity in the 1-d variational solution. It is 
shown that relatively accurate temperature and moisture retrievals can be achieved below 
optically thin clouds. For optically thick clouds, accurate temperature and moisture profiles 
down to cloud top level are obtained. For both optically thin and thick cloud situations, the 
cloud top height can be retrieved with relatively high accuracy (i.e., error < 1 km). NPOESS 
Airborne Sounder Testbed - Interferometer (NAST-I) retrievals from the Atlantic-THORPEX 
Regional Campaign are compared with coincident observations obtained from dropsondes and 
the nadir-pointing Cloud Physics Lidar (CPL). This work was motivated by the need to obtain 
solutions for atmospheric soundings from infrared radiances observed for every individual field 
of view, regardless of cloud cover, from future ultraspectral geostationary satellite sounding 
instruments, such as the Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) and 
the Hyperspectral Environmental Suite (HES). However, this retrieval approach can also be 
applied to the ultraspectral sounding instruments to fly on Polar satellites, such as the Infrared 
Atmospheric Sounding Interferometer (IASI) on the European MetOp satellite, the Cross-track 
Infrared Sounder (CrIS) on the NPOESS Preparatory Project and the following NPOESS series 
of satellites. 
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1. Introduction 


Nadir observations from a spacecraft- or an aircraft-flown infrared instrument can be 
used to infer the atmospheric temperature, moisture, and concentration of other chemical species 
using radiative transfer equation inversion techniques. The retrievals of atmospheric state, 
temperature and moisture profiles, obtained from infrared radiometric measurements will contain 
intolerable error near and below the cloud level if the cloud radiation and the attenuation of 
infrared radiation emitted from the Earth’s surface and the atmosphere below the clouds is not 
properly accounted for in the retrieval process. Clouds greatly complicate the interpretation of 
infrared sounding data. The new ultraspectral resolution infrared spectral radiance data contain 
the information needed to alleviate much of the ambiguity between cloud, atmospheric 
temperature, and moisture contributions which exists in lower spectral resolution sounding 
radiance data. Since there are vast cloudy regions of the globe, a great deal of effort has gone 
into cloud detection and cloud-clearing processes (Smith et al. 2004). Nevertheless, the schemes 
dealing with cloud detection and cloud-clearing (Smith 1968) remain a major source of error in 
the final retrieval products. Some schemes limit themselves to dealing with the observations 
unaffected by clouds (e.g., Chedin et al. 1985), while others make direct use of the cloudy 
radiances and attempt to retrieve temperature and moisture along with the cloud parameters (e.g., 
Susskind et al. 1984). Recently, fast molecular and cloud transmittance models (Huang et al. 
2004) have been developed and used to retrieve cloud optical thickness from Atmospheric 
InfraRed Sounder (AIRS) measurements (Wei et al. 2004). The retrieval methodology is reliant 
on other cloud and thermodynamic parameters, such as cloud top pressure and atmospheric 
profiles of temperature, moisture, and ozone. Fast molecular and cloud transmittance models 
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have also been used to enable the infrared radiances to be used under cloudy conditions with the 
accuracy required for sounding retrieval processing; the EOF (empirical orthogonal function) 
statistical regression retrieval algorithm (e.g., Smith and Woolf 1976; Zhou et al. 2002) has been 
expanded to include realistic cloud parameters (e.g., cloud top height, effective particle diameter, 
and optical thickness) to deal with cloudy as well as cloud-free observations (Smith et al. 2004; 
Zhou et al. 2005a). With that recently developed EOF regression algorithm, cloud parameters as 
well as atmospheric profiles are retrieved simultaneously from infrared (IR) spectral radiance 
observations. 

The NPOESS (National Polar-orbiting Operational Environmental Satellite System) 
Airborne Sounder Testbed - Interferometer (NAST-I) has been successfully operating on high 
altitude aircraft since 1998 (e.g., Cousins and Smith 1997; Smith et al. 2005). NAST-I is 
designed to support the development of future satellite temperature and moisture sounders such 
as the Infrared Atmospheric Sounding Interferometer (IASI) on the MetOp satellite, the Cross- 
track Infrared Sounder (CrIS) on the NPOESS Preparatory Project (NPP) and the following 
NPOESS series of satellites, as well as the GIFTS and the HES designed to fly on geostationary 
satellites. For obtaining dynamics (i.e., moisture flux and winds) from geostationary 
ultraspectral infrared radiance data, atmospheric profiles must be retrieved for every field of 
view, regardless of cloud cover, in order to obtain motion measurements from a time sequence of 
the 3-dimensional images of the atmospheric state variables (i.e., contiguous and continuous 
retrieval data are needed to construct the 3-dimensional images). Here we report that a physical 
retrieval scheme is developed to further improve retrieval accuracy based on EOF regressions. 
The physical retrieval scheme is described alone with the retrieval results for demonstration. 
Retrieval results of cloud and atmospheric properties from NAST-I observations are compared 
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with coincident observations obtained from the nadir-pointing Cloud Physics Lidar (CPL) and 
dropsondes, respectively. 

2. Physical retrieval scheme 

a. Radiative transfer model and Jacobian matrix 

The radiance measurements within the short wavelength region, where the observed 
radiance may be affected by reflected solar radiation, are typically not used during daytime 
observing conditions. The cloud transmissive and reflective functions can be coupled with a fast 
molecular radiative transfer model. The NAST-I fast transmittance model used here is a 
combination of the Optimal Spectral Sampling (OSS) fast molecular radiative transfer model 
(Moncet et al. 2003; Liu et al. 2003) and the physically-based cloud radiative transfer model 
based on the Discrete Ordinate Radiative Transfer (DISORT) (Stamnes et al. 1988) calculations 
performed for a wide variety of cloud microphysical properties (e.g., Yang et al. 2001; Huang et 
al. 2004). In the OSS approach, an extension of the Exponential Sum Fitting Transmittance 
method (Wiscombe and Evans 1977), radiance for each instrument channel is represented as a 
linear combination of radiances computed at a few pre-selected monochromatic frequencies 
within the domain spanned by the instrument line shape function (Moncet et al. 2003). Since the 
OSS model calculates channel radiances and transmittances at a few representative 
monochromatic frequencies, the top of atmosphere radiance for an atmosphere containing a 
single cloud layer is expressed by the monochromatic radiative transfer equation at those 
frequencies, 


R = R U F l X tc + R c X tc + + R l Fr X tc ■ 


( 1 ) 


5 



where R is the upwelling spectral radiance at the top of atmosphere. Ft and Fr are the cloud 
transmissive (including both direct and diffuse parts) and reflective (or albedo) functions, 


respectively. The terms R 0 , R c , R { , and R{ are upwelling emission below the cloud, emission 

from the cloud, upwelling emission above the cloud, and downwelling emission above the cloud, 
respectively. These terms are expressed below: 

i 

^0 = £ B s + JB dr + (1 - £ )RqT cs , (2) 

T c* 

< =r„ ^F T + i? c )f jBr/r', (3) 

i? f =(l-F R -F T )B(r c ), (4) 

1 

R\ = jBJr , ( 5 ) 

x tc 

1 

R } 1 =jBdr\ ( 6 ) 

r tc 

where e refers to the Earth’s surface emissivity, B is the Planck function, r is the total 
transmittance from any given level to an upper boundary such as cloud level or the top of the 
atmosphere, and r' is the total transmittance from any given level to a lower boundary such as 
cloud level or the Earth’s surface. The transmittance between the cloud level and the Earth’s 
surface is r cs , while r tc is the transmittance between the top of the atmosphere and cloud level. 
The subscripts s, c, and t denote the surface, cloud, and top of atmosphere, respectively. Figure 1 
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illustrates the various terms described above. It is noted that the top of atmosphere is assumed to 
be at the aircraft altitude for applications with the NAST airborne data. The transmissive and 
reflective functions at a given wavelength are taken from the pre-calculated database. In this 
study, particle habits are assumed to be hexagonal columns. Examples of cloudy radiance 
spectra simulated by the NAST-I forward radiative transfer model against cloud-free conditions 
are illustrated in Figure 2, showing different cloud features are captured mostly in the long 
wavelength window regions. 

The cloud phase (i.e., clear, liquid, and ice cloud) is simply defined according to the 
retrieval. In other words, retrieved cloud top temperature determines whether the cloud is liquid 
or ice. The relationship between the optical thickness, ice water path, and particle size is 
simplified by the relationship between the ice water path and cloud optical thickness based on 
the parameterization of the balloon and aircraft cloud microphysical database (Heymsfield et al. 
2003). Here, a simplified relation is used to specify cloud effective particle diameter from the 
cloud optical thickness, 


D = 


a£ w 

^ VIS 


(E . -b£ w ) 

' ~VIS ~VIS' 


(V) 


where a, b, and w are constants; D e is the effective cloud particle diameter; and Zg vis (later just g) 
is the visible optical thickness of the cloud. Two key cloud parameters, namely cloud top 
pressure and cloud optical thickness, remain in the regression and the first part of the physical 
iteration inversion. The weighting functions (or Jacobian matrixes) for cloud parameters are 
computed by a numerical perturbation method while others are computed by an analytical 
scheme. NAST-I analytical Jacobian matrixes for temperature, water vapor, and skin 
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temperature were computed as shown in Zhou et al. (2002). The numerical Jacobian matrix for 
cloud top pressure and optical thickness are expressed as 


A 


ij 


\R,(V’j)- R,(V’j 
H’j 


( 8 ) 


where i is channel index, / is cloud parameter ip index ( ipi is cloud top pressure and ip 2 is cloud 
optical thickness), and R is the simulated spectral radiance. An example of cloud parameter 
Jacobian is shown in Figure 3. The different spectral distributions of these weighting functions 
indicate the different channel radiance sensitivity and/or response to the cloud height and optical 
thickness. 


b. Inversion scheme 

An iterative eigenvector-regression-retrieval scheme was initially developed and used for 
retrieval analyses on both thermodynamic profiles and cloud parameters. These regression 
retrievals were compared with CPL and dropsonde data indicating a favorable agreement. Once 
the first guess is generated from the regression technique described by Zhou et al. (2005a), a 
non-linear iterative procedure is set up to produce a retrieval that is an improvement of the first 
guess (i.e., the EOF regression retrieval results). A one dimensional (1-d) variational solution, 
also known as the regularization algorithm or the minimum information method (e.g., Twomey 
1963; Tikhonov 1963; Rodgers 1976; Hansen 1998), is chosen for NAST-I physical retrieval 
methodology which uses the regression solution as the initial guess. This solution has been 
applied to NAST-I radiances for retrieving both thermodynamic parameters and trace gases 
under cloud-free conditions (Zhou et al. 2002; Zhou et al. 2005b). Here, we expand the same 
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methodology to include cloud parameters dealing with retrieval problems under both cloudy and 
cloud-free conditions. 

If the NAST-I observed radiance R™ of each channel is known, then R'" can be 

considered as a non-linear function of the atmospheric temperature profile ( T ), water vapor 
mixing ratio profile ( q ), surface skin temperature (T v ), surface emissivity (e), cloud top pressure 
(P c ), effective cloud particle diameter ( D e ), cloud visible optical-thickness (if), cloud phase (0), 
etc. Therefore R 1 " equals R^,(T,q,T s , a, P c ,D e £,(j),.. .) + a v , wherea v is the instrument plus other 
sources of noise. Notice that D e is initially a function of £ as in Eq. (7). In general, 


Y"'=Y(X) + cr, (9) 

where the state vector X contains atmospheric temperatures, atmospheric moisture mixing ratios, 
surface skin temperature, cloud optical thickness, cloud top height, etc. Y" ! contains N (number 
of channels used) observed radiances. The linear form of Eq. (9) is 

aY = Y'aX, (10) 

where Y' (or A) is the linear tangent of the forward model Y (or R) , the weighting function (or 
Jacobian) matrix. Here the linear model Y' uses an efficient analytical form (Li et al. 2000) for 
thermodynamic parameters and a numerical perturbation method for cloud parameters. A 
general form of the 1-d minimum variance solution minimizes the following penalty function 
(Rodgers 1976), 
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( 11 ) 


J(X) = [ Y m - Y(X)] t E' 1 [ Y“ - Y(X)] + [X - X 0 ] 1 H [X - X 0 ] , 


where superscript T denotes the transpose. Using the Newtonian iteration. 


x n+1 =x n+ r(x n )-‘j'(x n ), 


( 12 ) 


the following Quasi-Nonlinear iterative form (Eyre 1989) is obtained: 


aX n+1 =(Y: t E- 1 Y>H)- 1 K t E- 1 (aY n + KaX n ), 


(13) 


where aX n = X n - X 0 , aY n = Y m - Y(X n ) , X is the thermodynamic and cloud parameters to be 
retrieved, X 0 is the initial state of these parameters or the first guess, Y" ! is the vector of the 

observed radiances or brightness temperatures used in the retrieval process, E is the observation 
error covariance matrix which includes instrument noise and forward model error, and H is the a 
priori matrix which constrains the solution. Here, H can be the inverse of the a priori first guess 
error covariance matrix or another type of matrix. If the statistics of both the measurement and a 
priori error covariance matrix are Gaussian, then the maximum likelihood solution is obtained. 
However, if the a priori error covariance matrix is not known or is estimated incorrectly, the 
solution will be suboptimal (Eyre 1989). Here we apply H = dl in Eq. (13), where y is a 
Lagrangian multiplier which serves as a smoothing factor. Eq. (13) becomes 


aX n+I = (Y^E-'K + aiy l v; 1 E“‘ (aY n + y; aXJ, 


(14) 
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which is commonly referred to as the minimum information solution. It is noted that y is 
dependent upon the observations, the observation error, and the first guess of the atmospheric 
profile; often it is chosen empirically (e.g., Susskind et al. 1984; Smith et al. 1985; Hayden 
1988). In the NAST-I retrieval procedure, the Discrepancy Principle (e.g., Morozov 1984; 
Carfora et al. 1998; Li and Huang 1999) is applied to determine the appropriate smoothing factor 
y . Thus 


Y[X(y)]-Y'" 2 =cr 2 


(15) 


N 

where o 2 = ^ e 2 k , e k is the square root of the diagonal of E or the observation error of channel 

k , which includes instrument error and forward model error. Error is shown in the equation, 
e\ = r\l + / k 2 , where rj k is the instrument noise of channel k while / k is the forward model error 

which is assumed to be 0.3 K for the same channel. Usually o 2 can be estimated from the 
instrument noise and the validation of the atmospheric transmittance model used in the retrieval. 
Since Eq. (15) has a unique solution for y (Li and Huang 1999), the atmospheric parameters and 
the smoothing factor can be determined simultaneously. In NAST-I retrieval processing, a 
simple numerical approach is adopted for solving Eq. (15); y is changed in each iteration 
according to y n+1 = q n y n , here q is a factor for y increasing or decreasing. Based on Eq. (15), 
q is obtained in each iteration by satisfying the following conditions: 


<7o=LO; 


if 


Y(X n )-Y m < a 2 , then q n =\.5\ 
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if 

if 


Y(X n ) - Y” 
Y(X n ) - Y n 


>o 2 , then q n = 0.5; 

= o 2 , then stop the iteration. 


The q factor has been found from empirical experience to insure that the solution is stable 
between iterations. Thus, y keeps changing until the iteration stops. 

In the retrieval processing, several checks are made for retrieval quality control. The 
RMS (root-mean-square) of quantity [Y(X;) - Y m ] from all selected channels X\ is computed to 
check the convergence (or divergence). If Xi+i > X\ within 2 iterations (i.e., iteration diverges), 
then the iteration is stopped and the retrieval is set to the first guess (or the previous atmospheric 
state); otherwise, iteration continues until ^ n <1.0K and \x n ~ Xn-i | <0.0 IK, or a maximum of 10 
iterations is reached. The degree of convergence of each iteration depends on the accuracy of the 
previous atmospheric and surface state. In addition, at each iteration, each level of the water 
vapor profile is checked for super-saturation. A unity magnitude of relative humidity is assumed 
at any super-saturated level. 

The cloud particle size is derived using the relationship of Eq. (7) between the particle 
size and optical thickness. This relationship plays an influential role in both regression retrieval 
and physical matrix inversion. However, this relationship may not be true from case to case; on 
the other hand, the cloud signatures are captured in the long wavelength window region as shown 
in Figure 2, illustrating that the radiance spectral slope is sensitive to particle size and radiance 
magnitude is sensitive to optical thickness. Therefore, cloud microphysical parameters, namely 
effective particle diameter and visible optical thickness, are further refined with the radiances 
observed within the 10.4 pm to 12.5 pm window region, similar to what was described by Wei et 
al. (2004). The fitting of the magnitude and spectral slope of the radiance spectrum within the 
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window region is important for refining the cloud microphysical parameters. If these cloud 
parameters vary by more than 10% during this stage, the matrix inversion is performed again to 
produce atmospheric temperature and moisture profiles. To summarize this hybrid inversion 
scheme, a simplified flowchart is shown in Figure 4. 

3. Case study and validation 

NAST-I instrumentation, measurements, calibration, and radiance validation are 
documented elsewhere (e.g., Cousins and Smith 1997; Smith et al. 1999; Larar et al. 2002; Smith 
et al. 2005). NAST-I provides relatively high spectral resolution (0.25 cm" 1 ) measurements in 
the spectral region of 645-2700 cm" 1 with moderate spatial resolution (a linear resolution equal 
to 13% of the aircraft altitude at nadir) cross track scanning. While a large amount of data has 
been collected since 1998 under a variety of meteorological conditions, results from only a very 
limited data set are presented herein for the purpose of “cloudy” retrieval demonstration. 
Retrievals from the Atlantic-THORPEX Regional Campaign (ATReC) (e.g., Shapiro and Thorpe 
2004) are used to demonstrate this inversion methodology. These data, together with the 
radiosondes and dropsondes released from the NOAA G-4 aircraft that flew below the NASA 
ER-2 aircraft, provide a unique data set for detailed analysis of retrieval resolution and accuracy. 
During this field campaign, cloud properties were also provided by the nadir-pointing Cloud 
Physics LIDAR (CPL) on board the NASA ER-2 aircraft (McGill et al. 2002). All coincident 
observations obtained during this experiment are used to understand the atmospheric state and 
cloud microphysical properties for validating NAST-I retrievals. 

The experiment of 5 December 2003 is chosen to test and demonstrate this inversion 
scheme with a realistic cloud radiative transfer model. The target scenes (latitude from 32°N to 
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42°N, longitude from 68°W to 76°W), shown in Figure 5, covered a variety of conditions desired 
by the scientific objectives of the experiment. These included a variety of cloud conditions, such 
as medium-level altocumulus, as well as low-level cumulus, thunderstorm, and extensive high 
cirrus in the ATReC region covered by the ER-2 and G-4 aircraft. The EOF regression retrievals 
have already shown reasonable agreement with the dropsondes and CPL observations (Zhou et 
al. 2005a). Here, the improvements through the physical retrieval scheme described above are 
emphasized by the demonstration of both radiance fitting and retrieval parameter validation. 

a. Physical retrieval results and validation 

The physical retrieval results of cloud and thennodynamic parameters are shown in 
Figure 6. NAST-I retrieved cloud top height ( H c ) from the nadir observations against CPL 
measured cloud top heights of the top 2 layers; and Figure 6b shows the cloud optical thickness 
inferred from NAST-I measurements against that of the CPL 1064 nm channel measurements. It 
is noted that NAST-I spatial resolution (at the cloud height and at nadir) is 13% of the distance 
between the aircraft altitude and the cloud height (i.e., 1.56 km when the cloud height is at 8 km 
and the ER-2 is at 20 km), while the CPL horizontal resolution is about 0.2 km; furthermore, the 
NAST-I vertical resolution is about 1 km while the CPL vertical resolution is 0.03 km. Despite 
the differences of the instruments and of their spatial resolutions, the cloud top heights inferred 
from NAST-I compare very well with CPL measurements for the variety of cloud conditions 
observed. The measurement sensitivity and accuracy of cloud optical thickness inferred from the 
infrared measurement is expected to be much poorer than that measured by the CPL because of 
the spatial resolution differences between the two instruments. Even so, NAST-I cloud optical 
thickness retrievals compare favorably to CPL observations. 
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NAST-I retrieved temperature and relative humidity (RH) vertical cross sections are 
shown in Figure 6d and 6e, respectively. The areas whited out are under the clouds where the 
cloud optical thickness is larger than one. The variation of atmospheric conditions is captured 
very well by NAST-I retrievals, not only for the clear regions above optically thick clouds, but 
also for regions below optically thin clouds. These soundings are also validated by the 
dropsondes released from the G-4 aircraft. The dropsondes are used to reveal the retrieval 
sounding accuracy under cloudy conditions. Inter-comparisons between each dropsonde and 
NAST-I local mean (i.e., 5 X 5 single field of views) retrieval are presented in Figure 7. In 
general, the retrievals show a good agreement above the clouds; the sounding comparison 
continues to show a good agreement under the optically thin clouds to the second layer clouds as 
is indicated by the CPL observations, and retrievals are relatively accurate under the optically 
thin clouds. The erroneous retrievals below opaque clouds are shown. These errors, which were 
due to the limitation of infrared observations, were expected and should be discarded. As shown 
in Figures 6 and 7, retrievals of temperature and moisture above the clouds are not disturbed by 
the clouds below. Both temperature and relative humidity profiles show a continuity variation 
from time to time (i.e., from location to location) indicating that the atmospheric features are 
well captured by NAST-I observations and retrievals. 

b. Improved with physical inversion 

The regression retrievals with realistic cloudy radiance training have already shown the 
improvement over the clear and/or isothermal equivalent radiance training (Zhou et al. 2005a). 
However, the physical inversion scheme further improves the retrieval accuracy from the cloudy 
iteration-regression. Figure 8 shows the deviation of retrieved parameters from the regression to 
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the physical algorithms used. Several approaches are used herein to show the retrievals are 
evidently improved through physical inversion. The cloud particle size is produced as a function 
of optical thickness in order to reduce a number of retrieval parameters in the regression process, 
which results in an artificial correlation between these parameters as in Eq. (7). Plotted in Figure 
9 are the retrieval results from regression (gray dots) to physical (black dots), the physical 
outcomes still obey the prediction of Eq. (7) but with expected scattering. 

The cloud parameters retrieved from NAST-I observations plotted in Figure 6 indicate 
that a variety of cloud conditions were also observed by the GOES-8 visible and infrared 
imaging sensors (Figure 5). The radiance fitting statistics over this large diversity data set gives 
a clear reflection on the integrity of this retrieval algorithm. A few thousand NAST-I 
observations (all nadir of 5 December 2003 flight) are used to compare with their associated 
retrieval simulations in order to show retrieval accuracy. The standard deviation error (STDE) 
and mean bias in brightness temperature are plotted in Figure 10; the retrieval accuracy is 
improved significantly through the physical inversion. The relatively large error in the short 
wavelength region (i.e., wavenumber greater than 2200 cm" 1 ) is mainly due to NAST-I 
instrument noise, solar component not accounted for in the simulation, and the cloud model was 
cut off at a wavenumber of 2500 cm" 1 . Despite the uncertainty of the radiative transfer model 
(especially in the short wavelength region), the STDE from the physical retrieval is more like the 
instrument noise estimated from instrument calibration. This improvement evidently indicates 
that retrieved parameters, in general, are accurately retrieved through this physical inversion 
scheme. 

Finally, spectral radiance comparisons of individual scans have also been examined over 
each dropsonde location where the temperature and moisture profiles can be validated to ensure 
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that the deviations shown in Figure 8 are in favor of approaching the real atmospheric conditions 
through the physical inversion. Two samples, the first and 7 th dropsondes indicated in Figure 6d, 
are shown in Figures 1 1 and 12, respectively. Again, the radiance and retrieved parameters were 
compared to show the integrity of the algorithm. Panels (a) and (b) illustrate the radiance 
convergence to NAST-I observation, while (c), (d), and (e) show the temperature and moisture 
convergence to the dropsonde through the physical inversion. The retrieved cloud parameters 
are indicated in the figure captions. 

4. Concluding remarks 

Our previous regression results indicate some success in the ability to retrieve 
information below scattered and partially transparent cirrus clouds, or clouds with an effective 
optical thickness of less than one. The initial EOF regression has laid a first step in dealing with 
infrared sounding data under cloudy conditions, which is now significantly improved by the 
physical iteration inversion described in this study. Results achieved with airborne NAST-I 
observations show that accuracies close to those achieved in totally cloud-free conditions can be 
achieved down to cloud top levels. The accuracy of the profile retrieved below cloud top level is 
dependent upon the optical thickness and fractional coverage of the clouds. Retrieval accuracy 
of temperature and moisture profiles is greatly improved by the physical inversion as shown by 
dropsonde validation. Therefore, the radiances can be accurately simulated by using physically 
inverted results, which is important for direct assimilation into a forecast model. The 
thermodynamic profile information might be obtained through a combination of cloud clearing 
and direct retrieval from the clouded radiances using a realistic cloud radiative transfer model. 
The correct implementation still requires a considerable research development effort; however, 
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cloudy sky radiative transfer models should enable the extraction of profile information from 
cloud contaminated radiances suitable for numerical weather prediction application. 
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FIGURE CAPTIONS: 


Fig. 1 . Schematic illustrating the contributions to the simulated monochromatic radiance. 

Fig. 2. Simulated NAST-I spectra with cloud-free (black) and different cloud conditions of ice 
cloud top pressure at 250 hPa: the blue curve indicates particle size ( D e ) of 80 urn and optical 
thickness (£) of 0.5; the green curve indicates D e of 80 urn and £ of 0.1; the red curve indicates 
D e of 10 pm and ^of 0.1. 

Fig. 3. The normalized absolute weighting functions of cloud top height {H c at 6.8 km) and 
optical thickness. 

Fig. 4. A simple flow diagram summarizing the physical retrieval scheme employed for cloudy 
and/or cloud-free spectral radiance. 

Fig. 5. GOES-8 infrared image (at 18:30 UTC on 5 December 2003) showing a variety of 
cloudy conditions in the region covered by the NASA ER-2 and the NOAA G-4 aircraft. The 
ER-2 flight track is plotted over the GOES image. The dropsondes from the G-4 aircraft are also 
marked with asterisks (and a D#). 

Fig. 6. Panel (a) plots NAST-I physically retrieved cloud top height compared with the CPL 
measured cloud top heights of the top 2 layers (LI and L2). Panel (b) plots NAST-I retrieved 
visible cloud optical thickness (COT or §) compared with the CPL measurement. Panel (c) plots 
NAST-I retrieved cloud particle size. Panels (d) and (e) plot NAST-I physically retrieved 
temperature and relative humidity vertical cross sections, respectively. The areas wiped off are 
under the top layer clouds where the cloud visible optical thickness is larger than one and under 
the lower “opaque” cloud. The black vertical bars in panel (d) indicate dropsonde locations. 

Fig. 7. Inter-comparisons between NAST-I local mean retrievals (red curve is above the cloud 
while green curve is below the cloud) and dropsondes (blue curve). 
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Fig. 8. The retrieved parameter deviates from the regression to physical retrievals, as 
demonstrated in the above plots of NAST-I physical minus regression retrievals. 

Fig. 9. The correlation between cloud optical thickness and effected particle diameter: gray dots 
(formed a smooth curve) are from the regressions and black dots are produced from the physical 
retrievals. 

Fig. 10. Spectral data consistency showing the retrieval accuracy (for all nadir observations of 5 
December 2003 flight): (a) standard deviation error and (b) mean bias (simulated minus 
measured). The black and gray curves are produced from physical retrievals and regressions, 
respectively. 

Fig. 11. Comparison of spectral radiance and retrievals of temperature and water vapor 
corresponding to cirrus ice clouds near the first dropsonde. The cloud parameters changed from 
regression are the H c of 5.9 km from 6.8 km, £ of 1.37 from 0.69, and D e of 42.8 pm from 51.4 
pm. 

Fig. 12. Comparison of spectral radiances and retrievals of temperature and water vapor 
corresponding to cirrus ice clouds near the 7 th dropsonde. The cloud parameters changed from 
regression are the H c of 3.7 km from 5.4 km, of 1.81 from 0.41, and D e of 63.4 pm from 44.7 
pm. 
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Fig. 1 . Schematic illustrating the contributions to the simulated monochromatic radiance. 
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Fig. 2. Simulated NAST-I spectra with cloud- free (black) and different cloud conditions of ice 
cloud top pressure at 250 hPa: the blue curve indicates particle size ( D e ) of 80 urn and optical 
thickness (§) of 0.5; the green curve indicates D e of 80 urn and £ of 0.1; the red curve indicates 
D e of 10 pm and £fof 0.1. 
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Fig. 3. The normalized absolute weighting functions of cloud top height ( H c at 6.8 km) and 
optical thic kn ess. 



Fig. 4. A simple flow diagram summarizing the physical retrieval scheme employed for cloudy 
as well as cloud-free spectral radiance. 
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Fig. 5. GOES-8 infrared image (at 18:30 UTC on 5 December 2003) showing a variety of 
cloudy conditions in the region covered by the NASA ER-2 and the NOAA G-4 aircraft. The 
ER-2 flight track is plotted over the GOES image. The dropsondes from the G-4 aircraft are also 
marked with asterisks (and a D#). 
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Fig. 6. Panel (a) plots NAST-I physically retrieved cloud top height compared with the CPL 
measured cloud top heights of the top 2 layers (LI and L2). Panel (b) plots NAST-I retrieved 
visible cloud optical thickness (COT or §) compared with the CPL measurement. Panel (c) plots 
NAST-I retrieved cloud particle size. Panels (d) and (e) plot NAST-I physically retrieved 
temperature and relative humidity vertical cross sections, respectively. The areas wiped off are 
under the top layer clouds where the cloud visible optical thickness is larger than one and under 
the lower “opaque” cloud. The black vertical bars in panel (d) indicate dropsonde locations. 
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Fig. 7. Inter-comparisons between NAST-I local mean retrievals (red curve is above the cloud 
while green curve is below the cloud) and dropsondes (blue curve). 
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Fig. 8. The retrieved parameter deviates from the regression to physical retrievals, as 
demonstrated in the above plots of NAST-I physical minus regression retrievals. 
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Fig. 9. The correlation between cloud optical thickness and effected particle diameter: gray dots 
(formed a smooth curve) are from the regressions and black dots are produced from the physical 
retrievals. 
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Fig. 10. Spectral data consistency showing the retrieval accuracy (for all nadir observations of 5 
December 2003 flight): (a) standard deviation error and (b) mean bias (simulated minus 
measured). The black and gray curves are produced from physical retrievals and regressions, 
respectively. 
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Fig. 11. Comparison of spectral radiance and retrievals of temperature and water vapor 
corresponding to cirrus ice clouds near the first dropsonde. The cloud parameters changed from 
regression are the H c of 5.9 km from 6.8 km, £ of 1.37 from 0.69, and D e of 42.8 pm from 51.4 
pm. 
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Fig. 12. Comparison of spectral radiances and retrievals of temperature and water vapor 
corresponding to cirrus ice clouds near the 7 th dropsonde. The cloud parameters changed from 
regression are the H c of 3.7 km from 5.4 km, £ of 1.81 from 0.41, and D e of 63.4 pm from 44.7 
pm. 
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